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Abstract 

An excitation dynamics of new quantum systems of N equidistant energy levels in a mono- 
chromatic field has been investigated. To obtain exact analytical solutions of dynamic equa- 
tions an analytical method based on orthogonal functions of a real argument has been proposed. 
Using the orthogonal Legendre functions we have found an exact analytical expression for a 
population probability amplitude of the level n. Various initial conditions for the excitation 
of N- level quantum systems have been considered. 
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1 INTRODUCTION 

Theoretical investigations of the dynamics of the multiphoton vibrational excitation of polyatomic 
molecules by an IR laser radiation are based on the analysis of multilevel quantum systems as 
molecular models. Many other problems of nonlinear optics, spectroscopy and laser physics also 
lead to the consideration of the dynamics of multilevel systems. Analytical results obtained up to 
now for multilevel systems dynamics do not cover all the interesting and necessary cases. 

In this work an analytical method has been presented that offers new possibility to obtain exact 
solutions describing the excitation dynamics of a number of new multilevel quantum systems in an 
infrared laser field. The method uses orthogonal functions to obtain dynamic equations solutions. 
It generalizes an analytical approach that is based on integral transform and orthogonal polyno- 
mials. A particular example with the use of the orthogonal Legendre function is given. Analytical 
solutions for two different TV-level quantum systems with equidistant levels are found. 
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2 ANALYTICAL METHOD 

The dynamics of quantum systems in the monochromatic infrared radiation Ei cos(w;i) is described 
by the equations 



da n (t) 
dt 



/ n+ l e"^ 1 * a n+1 (t) + f n e 1 ^ 1 a n _i(i) , 

(1) 



a n (i = 0) = <5„, m , n = 0,1,2,... 

for the population probability amplitudes a n (t) of the level n. Here t = At is the dimensionless 
time ( A = ^to,i Ei/(2h) is the Rabi frequency, r is ordinary time). Eqs. (1) are obtained from 
the Schrodinger equation for a multilevel system with the use of the rotating wave approximation. 
In such description the multilevel system has two its characteristics: the dipole moment /in-i.n 
and the dimensionless frequency detuning e n for every transition n — 1 n . The dipole moment 
function /„ describes the dependence of the radiative transitions moments 

fn (2) 

on the transition number (on energy), fi i is the dipole moment of the lowest transition. The 
detuning 

£« = Kn-i - ^/)/A (3) 

characterizes the difference of the frequency of the transition n and the laser radiation frequency. 

Eqs. (1) for the iV-level quantum system can be solved by various methods. In work [Q 
analytical method has been developed on the basis orthonormal polynomial sequences p n {x)/d n 
(d n is a norm). The polynomials satisfy the recurrence formula 

+ f n = {rx + s n ) . (4 ) 

It is supposed that the coefficients /„ in formula (4) and in eq. (1) are identical functions, the 
coefficients s n and the frequency detuning are connected in the following way 

En = s n s n — 1 • (5) 

A number of exact analytical solutions with the help of the expression 

B 

a n (t) = [ cr(x) Pm j- X ^ Pn j x ^ exp[it(rx + s n )} dx (6) 



A 

has been obtained. 

Thus, the polynomials p n {x) are orthogonal on the interval {A, B) with respect to the weight 
function a(x) and give rise to multilevel quantum systems. The coefficients /„ , s n , r of 
recurrence formula (4) define all characteristics of this quantum systems, i.e. the dipole moment 
function /„ (2), the frequency detuning e n (5) and spacing of the energy levels 

E„ = E + nhwi + hA(s n - s ) ■ (7) 
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3 GENERALIZATION OF ANALYTICAL METHOD 



Let ip a (z) are some orthogonal functions of one variable z (in special case, of course, they can 
be orthogonal polynomials). In contrast to polynomials p n {x) which are orthogonal with the 
Kronecker delta <5 m ,„ , the functions ip a (z) can be orthogonal with both the Kronecker delta and 
the Dirac (5-function. The scalar product 

D a . p = J <p a {z) <p}(z) dM(z) = , a + p (8) 

L 

can have infinite values. Here L is an one-dimensional region of the integration, the measure 
dM(z) of the integral is non- negative, a and (3 are real numbers, not only integer ones. The norm 
d a = y D a jj =a of the orthogonal function can be infinite. 

In order to obtain the dynamic equation solution we use orthogonal functions (p n (z) which 
satisfy the recurrence formula 

Jn+i —j V Jn —j = [r 0{Z) + S n ) — . (9) 

"n+1 "n-l a n 

If 0(z) = x and (p n (z) are the polynomials p n (x) then (9) becomes (4). 
We shall show that any relation (r, A a , B a , C a are constants) 

r 9{z) ip a (z) = A a ip a+ i(z) + B a tp a (z) + C a ip a ~i(z) (10) 

for orthogonal functions tp a (z) can be reduced to form (9), when there exists a — ao for which 

C ao = . (11) 

Then (10) can be written as 

A ao f ao +i{z) = [r 0(z) - B ao ] if ao {z) , 

(12) 

(z) = [r 6{z) - B ao+n ] (z) + C ao+n (fi aa+ n-l(z) ■ 

Without any limitation we can take ao = . Now it is obviously that the function ip n (z) has the 
form 

<p n (z) = (fio(z) Pn[d(z)} (13) 

where p n [8(z)] is a polynomial of the argument 0(z). 

Any orthonormal polynomial p n (x) /d n satisfies the recurrence formula (4) Q . The functions 

ip n (z) = (f (z) p n [6(z)]/d n (14) 

also satisfy formula (4) when x = 8(z) . And we obtain (9). Systems of functions with infinite 
norm d n have to be considered specially. 

Thus, the solution of problem (1) is found in the following way: 
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(a) with the help of orthogonal functions ip a - ao (z) (8) we choose such a sequence 
of functions ip n (z) , a — oq = n = 0, 1, 2, . . . , which form the complete basis 
of solutions of the system of the relations (12); 

(b) the recurrence formulae (12) reduce to pattern (9); 

(c) the coefficients 



a n{ t) = / ^1 ^1 e «^*(*)+'») dM(z) 



(15) 



L 



of the expansion of the integral transform kernel 

CA( t)Z ) = ^^£) e *(rfl'W+-») (i 6 ) 

to the series of the functions (p n (z) are the solutions of problem (1). 

If the expansion of kernel (16) to series of orthogonal functions (polynomials) tp n (z) is known, the 
calculation of integral (15) can be avoided. 

4 THE ORTHOGONAL LEGENDRE FUNCTIONS 



Let's give a particular example of orthogonal functions which satisfy recurrence formula (9). It's 
known || that there are sy 
Legendre functions P^(z) 



known 0] that there are systems of functions tfm\z) among the complete family of the 1st kind 



V £\ Z ) = (-l)^-m- 1 T (N - riP&^z) 

orthonormalized on [—1, 1] with the Kronecker 5-function 



(17) 

fi—m 



1 

ip W {z ) tp \ N \ z) -^ = S m , l d 2 m , (18) 



-l 



m and I are non-negative integer number. N > 1 is an integer parameter (if m > N , the right part 
limit (17) is calculated at /i — * to). The norm of these functions 

d m = {{N + to — 1) ! (N — to — 1)! /to} 1/2 (19) 

is infinite when to — . Therefore the reducing known recurrence formulae Q] for the functions 
PJ?(z) to relation (9) requires the use of the linear dependence properties of the Legendre functions 
with the integer parameters fi = to and v = n . As a result, the functions (17) satisfy the relation 

frn+1 — T 1 h fm —j = r 0{z) — . (20) 
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Here d m is norm (18), 9{z) = z/^/z 2 — 1 , and 



(TV — to) (TV + to- 1)1 1/2 

, //t u , 

(21) 



^ TO (TO — 1) J 



= 2{(TV-2)(TV+l)/2}' 



'1/2 



As in (1), the coefficient f m at to = is equal zero. Since f m =N = , we have the system of TV 
linearly independent recurrence formulae of form (9) (at s n = ). In this work the relation (20) 
for values to > TV is not considered. 

5 TV-LEVEL QUANTUM SYSTEMS 

Let's use functions 

m>0, JV>1, (22) 

which form the complete system of solutions of TV equations (20) (when to < TV ), in order to 
obtain new solutions of problem (1). 

Let's integrate expression (15) written for functions (fim\z) ( m < TV ) which are orthogonal 



on [-1,1] if 6{z) = z/^z 2 ^! 



frn (z) <fi y n (z) irt (z/V?*=T) dz 



2 m-n+l y^f | m (_/V + TO — 1)! (TV — TO — 1)! } 2 X 



n(TV + n-l)!\ 2 irt {irt) n+m ^ 



(TV-n-1)! J n! 



r(TV - fc - |) (^irt)"-"- 1 -* 

(TV + to - k - 1)! (TV - to - fc - 1)! jfc! T(TV + n - k) 



2 F 2 ^N + n, n+^: TV + n-k , 2n + 1 ; -2irf J . (23) 

This solution describes the resonant excitation of a TV-level quantum system with the equidistant 
energy spectrum 

E„ = E + nhuJi (24) 

where E is the zero level energy. The dipole moment of this system depends on n according to 
(21). At the initial time moment t = only m-level is excited. 

Let's obtain another solution of problem (1) for another dependence on n of the dipole moment 
function. We can always consider reverse numeration of levels for finite level quantum systems, 
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because the choice of "upper" and "lower" level is arbitrary in such a case. Let's renumber the 
levels of the quantum system: n — > (N — 1 — n) . 

As a result, we obtain an analytical solution with the help of formulae (21) and (23) 

a n (t, N) = 2 m+n - N + 2 V5F {m(N + m-l)\ (N - m - 1)!}' x 



(N - n - 1) (2N - n - 2)! \ 2 Jrt (irt) 

JisT-n-iy. 



N— n+m— 1 
e"° ~"7TT TTr" X 



" 1 r(JV-A - ±) (-2irt) Ar -™- 1 - fe 



^ (AT + m - k - 1)! (N - m - k - 1)! fc! T(2N -n-k-l) 



2 F 2 (^2N - n - 1, N - n - ^; 2N - n - k - 1, 2N - 2n - 1; -2ir/j (25) 

which describes the resonant excitation of a TV-level quantum system also with the equidistant 
energy spectrum, but with different dipole moment function 



f ri "( 2jV -^+ 1 ) ] 1/2 n<N 



(TV - n + : 

(26) 

r = {(iV-l)/2}- 1/2 . 

Thus, new analytical method for obtaining exact solutions of the problem of the radiative 
excitation of multilevel quantum systems has been proposed. It allows to model the excitation of 
systems with more complex dynamics for the description of which both orthogonal polynomials 
and orthogonal functions have to be used. With the help of the orthogonal Legendre functions 
exact solutions for two various A^-level quantum systems have been obtained. 
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